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RESUMEN 

El metodo de los elementos finitos tiene una gran importancia en el 
campo de la ingenierfa, ya que permiten dar solucion a problemas que 
practicamente son imposibles de resolverse con metodos matematicos 
convencionales. Para lo cual este metodo aplica como base fundamental una 
discretization del elemento en un numero finito de sectores que posean una 
geometrfa mas simple, los cuales a su vez se encuentran interconectados a traves 
de nudos. 

Se presenta el marco teorico que conduce al calculo de la matriz de rigidez 
del elemento finito Q4 suavizado por flexion y se desarrollan programas que pasan 
a enriquecer el sistema de computation CEINCI-LAB. Se realiza el analisis 
estatico y dinamico de varias vigas en voladizo: de hormigon armado y de acero 
con perfil tipo “I" 


ABSTRACT 

The finite element method has a great importance in the engineering field 
because it provides solutions to problems that are virtually impossible to resolve 
with conventional mathematical methods. For which this method applies as the 
foundation element discretization in a finite number of sectors with a more simple 
geometry, of which are interconnected by nodes. 


The theoretical framework lead to the calculation of the stiffness matrix 
Finite Element Q4 softened by bending and programs that enrich pass CEINCI- 
LAB computer system developed is presented. Static and dynamic analysis of 
cantilever beams make: reinforced concrete and steel profile type " I" 
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1. INTRODUCTION 

El metodo de elementos finitos se utiliza como tal no hace muchos anos, 
pero su desarrollo y conceptos ya se empleaban hace ya varios siglos. Entre los 
vestigios mas antiguos que dan fe de este tipo de calculo son claramente las 
piramides egipcias, ya que sus constructores para determinar su volumen 
utilizaron metodos de discretizado. 

Este metodo consiste basicamente en la division de un elemento continuo 
en un conjunto de elementos mas pequenos conectados entre si mediante nodos. 
El sistema discreto a ser utilizado esta constituido por estas divisiones, las cuales 
generan un sistema de ecuaciones que nos proporcionaran el estado de tensiones 
y deformaciones. Las propiedades de los elementos que unen a los nodos las 
determinan la geometrfa y material asignado a cada uno de ellos, y esto define su 
rigidez basado en las leyes de la elastica. La relacion entre fuerzas externas y las 
deformaciones se da basicamente por las relaciones constitutivas del elemento y 
la rigidez. (Onate 2013, Aguiar y Mroginski, 2010, Chandrupatla etal. 1999, Young 
etal. 2000) 


2. DESCRIPCION DEL ELEMENTO FINITO Q4 

El elemento finito Q4 es un elemento cuadrilatero, el cual dispone de dos 
grados de libertad por nudo, siendo estos un desplazamiento horizontal y un 
desplazamiento vertical respectivamente, adicionalmente desplazamientos en el 
centra para poder considerar el comportamiento a flexion. 

El elemento se encuentra en un sistema coordenado cuyo origen esta en 
el centra del elemento (coordenadas naturales). Las ordenadas de la elastica se 
encuentran gobernadas por las siguientes expresiones que indican el 
desplazamiento de un punto cualquiera del elemento finito. 


u(s, t ) = ul01 + u202 + u303 + u404 + u505 + u606 
v(s, t ) = 17101 + v202 + v303 + v404 + 17505 + v606 



Figura 1 . Coordenadas reales y naturales de un elemento finito cuadrilatero 
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Figura 2. Ordenadas de la Elastica (Fuente: Aguiar, apuntes 2009) 


El elemento finito tiene dos sistemas de coordenadas, las reales cuyo origen 
coincide con el origen del sistema coordenado X,Y y las coordenadas naturales 
cuyo origen se encuentra en el centra finito, como se puede ver: 
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Figura 3. Coordenadas reales y naturales de un elemento. 


2.1 Funciones de Forma 


Funcion 0, 



Figura 4. Funcion de forma 0 ^. 
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Funcion 0? 
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02 (s) = 
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Funcion 0 4 
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Figura 7. Funcion de forma 0 4 . 
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1 + 1 

04 (t) * v4 = — — — * v4 


04(t) = 


1 + t 


04 (s,t) = 0.25* (l-s)(l + t) 


Funcion 0g 



Figura 8. Funcion de forma 0 5 . 
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Figura 9. Funcion de forma 0 6 . 
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2.2 Derivadas de las funciones de forma 


( 7 ) 


A continuation se muestran las funciones de forma que sirven para 
interpolar la geometna y el campo de desplazamientos. 

Tabla 1 Derivadas de las funciones de forma con respecto a “s”, “t 


Derivadas con respecto a s 

Derivadas con respecto a t 

90 1 (s, t) 

g s =-0.25(1 -t) 

90i(s,t) 

at = °' 25(1 s) 

d0 2 (s, t) 

8s =0-25(1-0 

90z(s,O , 

at = °- 25(1 + !) 

a03 as' t) = o - 25(1 + t) 

d0 3 (s,t) nnr „ , x 

at = 0 ' 25(1 + s: > 

904 (s, t) 

a = 0.25 (1 + 0 

OS 

904 (S,0 

at = 025 (1 

90s (s, 0 

, = 2s 

ds 

90s (s, 0 _ 0 
dt 

90 6 (s,0 _ Q 
ds 

906 (s, 0 

dt 


2.3 Relacion entre deformaciones y desplazamientos 

De la teorfa de la Elasticidad se sabe que las deformaciones unitarias s X: e y en 
sentido x, y; la distorsion angular y xy valen: 
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£ = 


du 
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dv 

£y ~dy 
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( 8 ) 


Al reemplazar las ecuaciones en las deformaciones y distorsion angular, se 
tiene luego de realizar las respectivas derivadas, lo siguiente: 
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2.4 Matriz B 


B es la matriz gradiente de deformacion del elemento y Bi es la submatriz 
asociada al nodo i. Tiene la siguiente forma: 


I - d0j(s,i) 

I dx 


Bt = 


d0j(s,t) 
L dy 


0 


d0j(s,t ) 
dy 

d0j(s,t ) 


(i = 1,2, 3, 4, 5, 6) 


dx J 


( 10 ) 


La primera columna de la matriz B t esta relacionada con el 
desplazamiento u, la segunda con el desplazamiento v. Por este motivo el 
subfndice varfa de 1 a 6. Los grados de libertad del elemento finito se agrupan en 
el vector p. De tal manera que: 

= [ul ul u2 u2 u3 u3 u4 u4 u5 u5 u6 u6] 

(ID 


Por lo tanto, la relacion entre las deformaciones s con las componentes de 
desplazamiento p en funcion de las matrices B t viene definido por 

s = [B 1 b 2 b 3 b 4 B 5 B 6 ] p 


£ = B p 


( 12 ) 
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3. MATRIZ JACOBIANA 


Denominada simplemente Jacobiano relaciona las derivadas cartesianas y 
naturales de una funcion. Aplicada a las funciones de forma se tiene: 


j = 


dx (s,t) dy (s,t) 


ds ds 

dx ( s,t ) dy (s,t) 


ds 


ds 


(13) 


dx (s, t) 


ds 


= 0.25[— (1 - t)x 1 + (1 - t)x2 + (1 4- t)x 3 - (1 - t)x 4] 


dy(s,t ) 


dt 


= 0.25[ — (1 - t)y 1 - (1 - t)y2 + (1 + t)y3 - (1 - t)y4] 


(14) 


dx(s, t) 

dt 

dy(s, t) 
dt 


0.25[— (1 - s)xl + (1 + s)x 2 + (1 + s)x 3 - (1 - s)x4] 


0.25[— (1 - s)y 1 - (1 + s)y2 + (1 + s)y3 - (1 - s)y4] 


Los puntos auxiliares corresponden a s=0; t=0. Por lo tanto al reemplazar 
estos valores en la (Ec.7) se tiene la matriz Jacobiana para los puntos auxiliares. 


dx (s, t) 

ds 

dy Q,t) 
ds 

dx ( s , t) 
ds 

dy ( 5 , t) 
ds 


0.25[— xl + x2 + x3 — x4\ 

= 0.25[— yl + y2 + y3 — y4] 

= 0.25[-xl - x2 + x3 + x4\ 
= 0.25[— yl — y2 + y3 + y4] 


(15) 
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30; 30; dx 30; dy 

3s dx ds + dy ds 

d</>i _ 30; dx 30; 3y 

dt dx dt 3y dt 


Escrito en forma matricial, se tiene: 


'dx dy' 
ds ds 
dx dy 
.dt dt. 
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dx 
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De donde la derivada de las funciones de forma respecto de las 
coordenadas reales puede obtenerse de: 


Sea la matriz / _1 de la forma: 
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5(1,1) = yi (1,1) a01 g^' t) +/ 1 (1,2) ^ = n (1,1) * [-0.25 * (1 - t(0)] +/I (1,2) * [-0.25 * (l - s(i))] 


5(2,1) = 0 
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5(3,1) =yi (2,1) ^ +/I (2,2) ^ = yi (2,1) * [-0.25 * (l - t(i))] +yi (2,2) * [-0.25 * (1 - s(i))] 


5(1,2) = 0 


5(2,2) =yi (2,1) (2,2)^ = yi (2,1) * [-0.25 * (l - t(i))] + /I(2,2) * [-0.25(1 - s(i))] 

6(3,2) =yi (1,1) ^ +yi (1,2) ^ =yi (i,i) * [-0.25 * (1 - t(o)] +yi (1,2) * [-0.25 * (1 - s(0)] 

Para las columnas 9 y 10, de la matriz de compatibilidad de deformaciones, se 
tiene: 

5(1,9) = JLAUX( 1,1) * [-2 * s(i)] 5(1,10) = 0 

5(2,9) = 0 5(2,10) = JLAUX( 2,1) * [-2 * s(i)] 

5 (3,9) = JLAUX(2,2) * [-2 * s(r)] 5(3,10) = JLAUX(1,1) * [-2 * s(i)] 

Finalmente, para las columnas 11 y 12, se tiene: 

5(1,11) = jLAUX (1,2) * [-2 * t( i)] 5(1,12) = 0 

5(2,11) = 0 5(2,12) = JLAUX( 2,2) * [-2 * t(i)] 

5 (3,11) = JLAUX(2,2) * [-2 * s(i)] 5(3,12) = JLAUX (1,2) * [-2 * t(i)] 

4. RELACION ESFUERZO DEFORMACION 

Matriz para Tension Plana usada cuando una de las dimensiones es 
menor que la longitud (espesor). 


C T = 


1 -v 2 


v 

1 


0 

0 


0 0 — 

2 . 


( 21 ) 


Donde E es el modulo de elasticidad del material, v es el modulo de 
Poisson. 


Al reemplazar los vectores de esfuerzo 0 (tension) y en vector de 
deformaciones se tiene: 


" ®x " 


'1 v 0 ' 


^x " 

Gy 

E 

1—v 2 

v 1 0 

_ ^ 1—v 


a y 

Jxy. 


0 0 — 

2 J 


-Vxy . 


( 22 ) 


La matriz de rigidez K, del elemento finito se halla a partir de la energfa 
potencial de deformacion. Para ello se escribe la energfa potencial en funcion de 
los desplazamientos que es igual a la energfa potencial de deformacion del 
elemento, esto es: 

p* K e p = J v £ t C T edV 


(23) 
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Donde dV es el diferencial de volumen. Al reemplazar la ecuacion (Ec.6) y 
luego de simplificar terminos se halla: 

Ke = 


Al efectuar la integracion en coordenadas naturales (normalizadas) se tiene: 

K e = L\ L\ B l C T Be dA (25) 

Donde e es el espesor del elemento que se considera constante y dA es el 
diferencial de area. Para la integracion, es suficiente considerar 4 puntos de la 
cuadratura de Gauss. 

Con lo que se tiene: 


Jv 


B 1 C T B dV 


(24) 


K e = Vi^B\C T Bi det(/ ) e P L (26) 


5. SUBRUTINAS DE CEINCI-LAB PARA EL ANALISIS 
ESTATICO CON ELEMENTOS FINITOS. 

• [CG, ngl, ne]=cg_ef_Q4(divx, divy): 

CG Matriz de coordenadas generalizadas 
ngl Numero de grados de libertad 

ne Numero de elementos finitos para una viga empotrada en voladizo 
divx Nro. Divisiones en la longitud de la viga 
divy Nro. Divisiones en la altura de la viga 

• [MC] =mc_ef_Q4(divx,divy): 

MC Matriz de conectividad para una viga empotrada en voladizo. 
divx Nro. Divisiones en la longitud de la viga 
divy Nro. Divisiones en la altura de la viga 

• [VC]=vc_ef_Q4(CG,MC,ne): 

VC Matriz de vector de colocacion de una viga empotrada en voladizo. 
CG Matriz de coordenadas generalizadas 
MC Matriz de conectividad para una viga empotrada en voladizo. 
ne Numero de elementos finitos para una viga empotrada en voladizo 

• [KT] = KT_ef_Q4(k3,VC,ne,ngl): 

KT Matriz de rigidez total de la viga empotrada en voladizo KT. 
k3 Matriz que contiene la rigidez de miembro 

VC Matriz de vector de colocacion de una viga empotrada en voladizo. 

ne Numero de elementos finitos para una viga empotrada en voladizo 

ngl Numero de grados de libertad 
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[XE,YE] =gn_ef_Q4(1_,h,divx,divy): 

XE, YE Coordenadas reales del elemento finito Q4 

L Longitud de la viga 

H Altura de la viga 

divx Nro. Divisiones en la longitud de la viga 

divy Nro. Divisiones en la altura de la viga 


• [k3]=rigidez_elemento_finito{espesor,XE,YE,E,po\ssor\): 

k3 Matriz que contiene la rigidez de miembro 
espesor Espesor de la viga empotrada en voladizo 
poisson Modulo de Poisson 
E Modulo de elasticidad 

XE, YE Vectores que contienen las coordenadas del elemento finito 

• [Q] = Q_e/LQ4(Fn,ngl,CG): 

Q Matriz de cargas de la viga empotrada en voladizo 

Fn Matriz que contiene el o los nudos de los elementos finitos 

cargados y sus cargas verticales y horizontales. 
ngl Numero de grados de libertad 

CG Matriz de coordenadas generalizadas 


5.1 Analisis estatico de una viga rectangular en voladizo sobre la 
que gravita una carga puntual. 


Datos para el calculo 

Longitud= 3.00 m 

Altura= 0.40 m 

Espesor= 0.30 m 

Numero de divisiones en “X”= 4 

Numero de divisiones en “Y”= 2 



Figura 10. Geometrfa de la Viga 
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Ubicacion de Carga 
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Figura 11 . Ubicacion de la carga. 
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Figura 12 . Numeracion de elementos. 

Numeracion de Nudos 



Figura 13 . Numeracion de nudos 

Numeracion de grados de libertad 
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5.2 Solucion 


Grafica desplazamiento vs numero de divisiones 



Grafico 1. Nro. De divisiones vs desplazamientos (Viga rectangular) 

Kv= 431.2253 (CEINCI-LAB) 


5.3 Analisis estatico de una viga metalica tipo “I” en voladizo 
sobre la que gravita una carga puntual. 

Para el siguiente ejemplo se tiene una viga metalica tipo “I” la cual ha sido 
dividida en veintiocho elementos finitos con una carga puntual de 2 T/m, la viga 
tiene una longitud de 3 m y su seccion se muestra a continuation. 



Datos para el calculo 

Longitud= 3.00 m 
Altura= 0.25 m 

Numero de divisiones en el alma: 
. “X”= 4 

. “Y”= 3 

Numero de divisiones en el ala: 
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. “X”= 4 

• “Y”= 2 


Seccion de la Viga: 



Figura 16. Grafico que contiene Secciones de la Viga 
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Figura 17. Grafico que contiene dimensiones de la Viga 
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Figura 18. Numeracion de elementos 
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Numeracion de grados de libertad 
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Figura 19. Numeracion de grados de libertad 


5.4 Solucion 

Grafica desplazamientos vs numero de divisiones 
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Grafico 2. Nro. De divisiones vs desplazamientos (Viga I) 

Kv= 46.3223 (CEINCI-LAB) 

6. SUBRUTINAS DE CEINCI-LAB PARA EL ANALISIS 
DINAMICO CON ELEMENTOS FINITOS. 

En este apartado, se presentan los programas que sirven para el analisis 
dinamico de estructuras y que se encuentran en Aguiar (2012). Se los indica para 
tener completo el marco teorico del analisis sismico utilizando el Metodo Espectral, 

• [T,phi,OM]=orden_e/g(K,M): 

KL Matriz de rigidez lateral 

M Matriz de masas. 

T Periodo de vibracion. 

Fi Modos de vibracion. 

OM Frecuencia 
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• [A6]=espectro_nec15_mod_DE(R,i\p,i\e,B,to,J) 

R Factor de reduccion de las fuerzas sismicas para superestructura 
tip Factor de irregularidades en planta 
fie Factor de irregularidades en elevacion 

B Factor de reduccion del espectro en funcion del amortiguamiento 

to Periodo a partir del cual se realiza la reduccion para aisladores 

T Periodos de vibracion 

Ad Aceleraciones Espectrales para cada periodo de vibracion en m/seg 2 

• [qt]=c/esp/azam/enfos_moc/a/es_CQC(T,phi,Ad,FP,na,Wn,z) 

T Son los periodos de vibracion en cada modo 

phi Es una matriz que contiene a los modos de vibracion 

Ad Es un vector que contiene las aceleraciones espectrales en cada 
modo de vibracion. 

FP Es un vector columna con los factores de participacion de cada 
modo 

na Numero de modos de vibracion o grados de libertad 
Wn Vector que contiene las frecuencias de vibracion 
z Factor de amortiguamiento 

q Es la matriz que contiene los desplazamientos en cada modo 
qt Es el vector de desplazamientos resultantes luego de la combinacion 
modal 

• [Vt,Ft]=ft/erzas_mo£/a/es_CQC(M,phi,Ad,gama,na,Wn,z 

M Es la matriz de masas 

phi Es una matriz que contiene a los modos de vibracion, 

Ad Es un vector que contiene las aceleraciones espectrales en cada 
modo de vibracion. 

gama Es un vector con los factores de participacion de cada modo 
na Numero de modos de vibracion. 

Wn Vector que contiene las frecuencias de vibracion 
z Factor de amortiguamiento 

P Es la matriz que contiene las fuerzas en cada modo. 

V Es la matriz que contiene los cortantes en cada modo. 

Vt Es el vector de cortantes resultantes luego de aplicar la 
combinacion CQC 

Ft Es el vector de fuerzas resultantes obtenido con los cortantes 
resultantes. 

6.1 Analisis dinamico de una viga rectangular en voladizo sobre 
la que gravita una carga distribuida. 

Para el siguiente ejemplo se tiene una viga la cual ha sido dividida en doce 
elementos finitos con una carga distribuida de 2 T/m, la viga tiene una seccion de 
30/40 cm, 4 m de longitud. Con la carga distribuida se obtiene la matriz de masas. 
Se desea encontrar las respuestas maximas probables empleando el Metodo 
Espectral si la estructura se halla en Quito, en un suelo, cuyos factores de sitio 
son: F a = 1.2; F d = 1.3; F s = 1.3, que corresponden a un suelo tipo “C” (NEC-15). 
Se considera un factor de reduccion de las fuerzas sismicas R = 2 
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Datos para el calculo 

Longitud= 4.00 m 

Altura= 0.40 m 

Espesor= 0.30 m 

Numero de divisiones en “X”= 4 

Numero de divisiones en “Y”= 3 









I 

o 

d 










1 i 1 ^ i 1 ^ i 1 ^ i 



1 1 1 1 T 1 1 1 1 T | | | | T | | | | r 
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Figura 20. Geometrla de la Viga. 


Carga distribuida para hallar la matriz de masas 
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Figura 21. Carga uniforme distribuida 


Con la carga distribuida se encuentra las cargas puntuales en las 
coordenadas principales y se halla la matriz de masas que en el presente caso es 
de 4 por4. 

Numeracion de Elementos 
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Figura 22. Numeracion de elementos. 

Numeracion de Nudos 
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Figura 23. Numeracion de nudos 


Numeracion de grados de libertad 
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Figura 24. Numeracion de grados de libertad 


Discretizacion de las Masas 



6.2 Solucion utilizando CEINCI-LAB 


Los resultados que se encuentran se presentan a continuation 
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Figura 26. Grafico de CEINCI-LAB 


Tabla No. 2 Coeficientes para el espectro de analisis sismico. (NEC-15) 


Datos para el analisis sismico 

Factores 

Valor 

Zona Sismica 

Sierra 

Tipo de Suelo 

C 

Valor de R 

2 

Factor "Z" 

0.4 

Factor de Sitio "Fd" 

1.3 

Factor de sitio "Fa" 

1.2 

Factor de comportamiento 

1.3 

inelastico del suelo “Fs” 



Figura 27. Espectro de Diseno Normativa NEC-15 
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• Propiedades dinamicas 


Tabla No. 3 Propiec 

lades dinamicas de la Viga Rectangular. 


Modo 1 

Modo 2 

Modo 3 

Modo 4 

Peri'odo (s) 

0.2547 

0.038 

0.0122 

0.0055 

Modos de Vibracion 

0.1617 

0.7315 

1.4018 

1 .5407 

0.5822 

1.5359 

0.6154 

-1.3503 

1.1542 

0.8635 

-1.4725 

0.8087 

1 .7896 

-1.1226 

0.6228 

-0.2215 

Factor de Participacion 
modal 

0.7526 

0.4098 

0.2383 

0.1587 

Aceleracion Espectral (m/S2) 

3.8886 

2.427 

1.844 

1 .6929 

Desplazamiento Lateral (cm) 

0.16 

0.56 

1.11 

1.72 


6.3 Calculo de una viga metalica tipo “I” en voladizo sobre la que 
gravita una carga distribuida. 

Para el siguiente ejemplo, se tiene una viga metalica tipo “I” la cual ha sido 
dividida en dieciseis elementos finitos con una carga distribuida de 2 T/m, la viga 
tiene una longitud de 4 m y su seccion se muestra a continuacion. Se desea 
encontrar las respuestas maximas de desplazamiento ante la misma accion 
sismica del ejemplo anterior. 

Datos para el calculo 

Longitud= 4.00 m 
Altura= 0.55 m 

Numero de divisiones en el alma: 

. “X”= 4 

• “Y”= 2 


Numero de divisiones en el ala: 

• “X”= 4 

• “Y”= 1 

Seccion de la Viga: 



o 

m 




d 


Figura 28. Grafico que contiene dimensiones de la Viga 
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Figura 29. Geometria de la Viga. 

Carga Uniforme considerada 


Figura 30. Ubicacion de la carga. 

Numeracion de Elementos 
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Figura 31. Numeracion de elementos. 
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Figura 32. Numeracion de nudos. 
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Numeracion de grados de libertad 
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Figura 33. Numeracion de grados de libertad. 


Discretizacion de Masas 


La discretizacion de masas es similar que para la viga rectangular. Ver 
figura 25. 


Datos para definir el espectro de diseno del NEC-15 

Para definir el espectro de diseno del NEC-15 se utiliza los mismos datos 
mencionados en la Tabla 2. 

Espectro de Diseno Normativa NEC-15 

Se presenta el grafico del espectro de Diseno. Ver figura 27. 

❖ Propiedades dinamicas de la viga. 

Tabla No. 3 Propiedades dinamicas de la Viga Rectangular 



Modo 1 

Modo 2 

Modo 3 

Modo 4 

Periodo (s) 

0.1054 

0.0201 

0.0091 

0.0064 

Modos de Vibracion 

0.1892 

0.9629 

1.581 

1.1989 

0.6134 

1.5208 

0.124 

-1.4818 

1.168 

0.691 

-1.3775 

1 .0774 

1 .7674 

-1.0875 

0.6981 

-0.3261 

Factor de Participacion 
modal 

0.7629 

0.4259 

0.2093 

0.0956 

Aceleracion Espectral (m/S2) 

4.3207 

2.2458 

1.9718 

1.902 

Desplazamiento Lateral (cm) 

0.03 

0.1 

0.2 

0.3 


7. CONCLUSIONES 

• El metodo de elementos finitos es una herramienta potente que nos 
permite trabajar con diferentes geometrfas y caracterfsticas de los 
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materiales, obteniendo de una forma precisa los estados de 
deformaciones y tensiones 

• Al utilizar el metodo de elementos finitos podemos llegar a la 
conclusion que al realizar la discretizacion se obtienen resultados mas 
precisos realizando la mayor cantidad de divisiones, esto se refleja en 
que existe cambios en los valores a partir del quinto decimal. 

• En el analisis de vigas en I es aconsejable el trabajar con divisiones 
diferentes en el sentido (y) para de esta manera hacer coincidir los 
grados de libertad que comparten en la union de las alas con el alma, 
es decir teniendo una geometrfa diferente en los elementos finitos de 
las alas con los elementos finitos del alma. 
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